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ABSTRACT 

We report on detections of the high-excitation CO J=6-5, J=4-3 hues in 
Mrk231, a prototypical Ultra Luminous Infrared Galaxy (ULIRG) and Seyfert 1 
QSO. These observations are combined with CO J=3-2, HCN J=4-3 (this work), 
and CO J=2-l, J=l-0, ^^CO J=2-l, HCN J=l-0 measurements taken from the 
literature to provide better constraints on the properties of the molecular gas in 
an extreme starburst/QSO in the local Universe. We find that the CO J=4-3 
and J=6-5 transitions trace a different gas phase from that dominating the lower 
three CO transitions, with n(H2) ~ (1 - 3) x 10^ cm'^ and Tk ~ (40 - 70) K. 
This phase is responsible for the luminous HCN emission, and contains most 
of the H2 gas mass of this galaxy. The total CO line cooling emanating from 
this dense phase is found similar to that of the [C ll] line at 158 /xm, suggest- 
ing a very different thermal balance to that seen in lower IR-luminosity galaxies, 
and one likely dominated by dense photon-dominated regions. Our dense "sam- 
pling" of the CO rotational ladder and the HCN lines enables us to produce well- 
constrained Spectral Line Energy Distributions (SLEDs) for the dense molecular 
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gas in Mrk231 and compare them to those of high redshift starbursts, many of 
which have SLEDs that may be affected by strong lensing. Finally, we use our 
local molecular line excitation template to assess the capabilities of future cm 
and mm/sub- mm arrays in detecting CO and HCN transitions in similar systems 
throughout the local and distant Universe. 

Subject headings: galaxies: individual (Mrk 231) — galaxies: ISM — galaxies: 
starburst — ISM: molecules 



1. Introduction 

Luminous and Ultraluminous Infrared Galaxies ((U)LIRGs) represent a truly remarkable 
class of objects. With far-infrared luminosities of Lpm > 10^^~^^Lq, these galaxies emit a 
significant fraction of their total bolometric luminosity at infrared wavelengths (Soifer et 
al. 1986; Sanders 1989). Although these rapidly evolving and often merging objects are 
frequently hosts of both powerful starbursts and AGN (Sanders & Mirabel 1996), there 
is increasing observational evidence to suggest that most of this IR and FIR emission is 
produced by extreme star formation events (Genzel et al. 1998, Downes & Solomon 1998). 
Observations suggest that the evolutionary track followed by ULIRGs once their molecular 
gas supply is exhausted will likely be towards field elliptical/lenticular galaxies of moderate 
mass (Genzel et al. 2001; Tacconi et al. 2002), rather than towards hosts of optically bright 
QSOs as originaly suggested by Sanders et al. (1988). 

Several studies of the molecular gas have been made since its ubiquitous presence in 
(U)LIRGs was first established (Tinney et al. 1990; Sanders, Scoville & Soifer 1991; Solomon 
et al. 1997; sec Sanders 1997 for a review). Intense far-UV radiation and strong tidal fields 
in these merging/starburst objects are likely to be responsible for the differing physical 
properties of the molecular gas relative to those prevailing in lower-power starbursts (Aalto 
et al. 1991, 1995; Casoh, Dupraz & Combes 1991), while interferometer maps reveal compact 
(~ (0.5 — 1) kpc) CO-emitting regions (Bryant 1997; Downes & Solomon 1998; Bryant & 
Scoville 1998). A large fraction of the total molecular gas mass in these objects is in a 
warm, dense gas phase (Solomon, Radford, & Downes 1990; Solomon, Downes, & Radford 
1992a), and the same type of gas may be also responsible for the bright high-excitation CO 
lines detected in several dust-enshrouded and more luminous starbursts and QSOs (LpiR ~ 
10^^ Lq) at high redshifts over the last decade (see Solomon & Vanden Bout 2005 for a 
review) . 

Detections of high-excitation CO transitions such as J=4-3 and J=6-5 in local galaxies 
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were first made in the early 1990s (e.g., Harris et al. 1991, Wild et al. 1992, Giisten et 
al. 1993). Poor atmospheric transmission under all but the very best observing conditions 
has prevented systematic surveys of these transitions in galaxies other than in local, compact 
(~ 500 pc) and low luminosity (LpiR ~ 4 x 10^° Lq) starbursts such as M82 (White et 
al. 1994; Mao et al. 2000, Wild et al. 1993, Ward et al. 2003), and the centers of a 
few nearby spirals (Nieten et al. 1999). The lack of large-format array receivers at these 
frequencies makes imaging of high excitation CO transitions difficult except with special 
instrumentation (Fixsen, Bennett, & Mather 1999; Kim et al. 2002). This is unfortunate 
since, when combined with lower-J CO and ^^CO transitions, CO J=4-3, J=6-5 line emission 
with E4/kB ~ 55 K, Eg/ke ~ 116 K and critical densities of n43 ~ 1.9 x lO^cm"^ and 
uqs ~ 6 X lO^cm"'^, can provide excellent diagnostics of the excitation conditions of the 
molecular gas. 

To fill this gap, we initiated a multi-transition CO, ^^CO and HCN survey of ~ 30 
LIRGs with which to undertake a detailed study of their molecular gas component, with 
particular emphasis on the warm and dense phase fueling their prodigious star formation 
rates. In this paper we report on the principal result of the ULIRG/QSO Mrk231, the 
first galaxy in our sample for which multi-transition observations have been completed. We 
present an assesssment of the physical conditions of its molecular gas reservoir based on 
an LVG analysis, which we then use to compute the cooling budget of the reservoir and to 
construct CO and HCN Spectral Line Energy Distributions (hereafter SLEDs). We compare 
the Mrk231 SLEDs to the typically poorly constrained SLEDs of high-z starbursts, and 
use them to assess the capabilities of both future cm/mm interferometers and ESA's high 
frequency spaceborne observatory Herschel in detecting these important line diagnostics in 
starbursts in the distant Universe. Throughout this work we assume a fiat A-dominated 
cosmology with Hq = 71kms~^ Mpc~^ and = 0.27. 



2. Observations and results 

The 15-meter James Clerk Maxwell Telescope (JCMT)E| on Mauna Kea in Hawaii was 
used for observations of the CO J=6-5, CO J=4-3, CO J=3-2 and HCN J=4-3 transitions in 
the archetypal ULIRG/QSO Mrk231. These observations are part of our large and ongoing 
CO, ^^CO and HCN line survey of 30 LIRGs (Lir > 10^^ Lq) which is being undertaken using 



^The James Clerk Maxwell telescope is operated by the Joint Astronomy Centre on behalf of the Sci- 
ence and Technology Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific 
Research and the National Research Council of Canada. 
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the JCMT and the IRAM 30-m telescope (see Papadopoulos et al. 2007 for a description). 

2.1. The CO J=6-5, 4-3 observations 

Receiver W at D-band (620-710 GHz) was used on February 20, 2005 in SSB mode to 
observe the CO J=6-5 hne ( z/j-est = 691.473 GHz) under excellent atmospheric conditions 
(^225 ^ 0.035), with typical system temperatures of Tgys ~ (4500 — 5200) K (including 
atmospheric absorption). The Digital Autocorrelation Spectrometer (DAS) was deployed 
in wideband-mode (1.8GHz), with a resulting instantaneous bandwidth of ~ 780kms~^ 
which was more than sufficient to cover the full width of the CO line (FWZI ~ 400kms~^). 
Rapid beam switching with a chop frequency of 1-2 Hz and a beam throw of 60" (in Az) 
yielded flat baselines after a total of 2 hours of integration (on-l-off). Pointing at D-band 
is complicated by the paucity of suitably bright, compact pointing sources. Differential 
pointing was therefore employed using a combination of B-band (330-360 GHz, B3 receiver) 
and D-band measurements, with the relative pointing offsets of the two different receivers 
determined prior to the observing run. Initial pointing and focus checks were made using 
Jupiter which, at the time of the observations, had a semi-diameter of ~ 20" - we note that 
pointing using sources that are extended relative to the primary beam has been shown to 
work well at these high frequencies. Pointing checks with B3 were made prior to, and during 
our observations, and were applied to the telescope pointing model for D-band, resulting in a 
typical rms scatter in the latter of ~ 2" — 3". The CO J=4-3 line (z/rcst = 461.0407 GHz) was 
observed with receiver W tuned to SSB mode in C-band (430-510 GHz) on April 23 2005. 
Dry conditions {T225 ~ 0.035) yielded typical system temperatures of Tsys ~ 1900 K. A DAS 
bandwidth of 920 MHz (~ 600 km s^^) was used, which was sufficient to cover the line with 
ample margin for baseline subtraction. Beam switching at 1-2 Hz with a beam throw of 
30" (Az) produced excellent flat baselines after a total of 40 mins (on+off) of observations. 
Pointing checks were again differential, with an rms uncertainty of ~ 3". 

Estimates of the aperture efficiencies at both frequencies were made from repeated 
observations of Mars, resulting in a value of ?7a(461 GHz) = 0.36 (for 6'hpbw = H") which is 
consistent within ~ 10% with values reported by others for the same observing periodj^ The 
scarcity of rja, measurements at 691 GHz did not allow such comparisons and so we adopted 
the mean value, r7a(691 GHz) = 0.18, of our two measurements (^hpbw ~ 9" from a Mars 
beam map). This value has a larger uncertainty - of around 30% - due to larger calibration 
uncertainties as well as the effects of thermal distortions/mechanical deformations of the dish 



^http: / /www .jach.hawaii.edu/JCMT/spectralJine/Standards/efF_web. html 
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affecting high-frequency observations even on prime sub-mm telescopes such as the JCMT. 

2.2. The CO J=3-2 and HCN J=4-3 observations 

The CO J=3-2 (345.7960 GHz) and HCN J=4-3 (354.5054 GHz) observations were 
made on July 10th 1999 (CO) and January 26th 2005, January 17th 2006 (HCN), with 
receiver B3 tuned in SSB mode with effective system temperatures of Tgyg = 670 K and 
Tgys = (390 — 470) K respectively. A DAS configuration with ~ 1.8 GHz bandwidth was used 
for CO J=3-2 (~ 1556 km s~^), whilst a narrower-band 920 MHz/dual-channcl mode was 
used for HCN J=4-3 (~ 776kms~^) for increased sensitivity. Beam switching at 1 Hz and a 
chop throw of 60" (Az) was employed in both cases, giving flat baselines for total integration 
times of 20min (CO) and 4 hours (HCN). A number of aperture efficiency measurements 
were made using Mars, yielding (r/a) ~ 0.515 ± 0.08 (for 6'hpbw = 14"). Frequent pointing 
checks were made by observing strong sources in both continuum and spectral line mode, 
yielding an rms uncertainty of ~ 3" for the pointing model residuals. Finally, observations 
of strong spectral hue standards such as IRC 10216, OMCl, W75N and W3(0H) were used 
to verify the amplitude calibration, to estimate calibration uncertainties (~ 15%), and to 
monitor the overall performance of the telescope. 

2.3. Data reduction, results 

All data were reduced using the JCMT spectral hue reduction package SPECX. Zero- 
order baselines were removed and spectra inspected individually, prior to being co-added to 
produce the final spectra shown in Figs. 1 and 2, where an excellent line profile agreement, 
with a FWZl ~ 400 km s^^ , is evident for all the transitions observed. The velocity-integrated 
line fiux densities were estimated from these spectra using 

Sune=/ S.dV = -^Ke(x) / TldV = l^:^M^^K.(x) / TldV, (1) 
Jay VaT^^ Jay % Jav 

where the term Kc(x) = x^/(l — e~^^), with x = ^s/(1-2^hpbw) and ^^^source diameter, 
accounts for the geometric coupling of the gaussian component of the beam with a finite- 
sized, disk-like source. For Mrk231 we have used a value of ~ 3" (likely to be even 
smaller for the high- J CO transitions and the HCN emission), obtained from interferometric 
observations of CO J=l-0 (Downes & Solomon 1998), yielding an upper limit of Kc(x) ~ 1.04 
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in the case of the CO J=6-5 hne. The hne fluxes obtained from observations presented 
in this paper, those extracted from observations reported in the hterature and brightness 
temperature hne ratios are hsted in Tables 1 and 2. 

2.3.1. CO J=l-0 line flux: single dish versus interferometers 

It can be seen from Table 1 that, with the exception of a single measurement by Kriigel 
et al., all CO J=l-0 single dish measurements are in good agreement with oneanother, 
but that they are higher than those obtained using interferometers (which are, in turn, 
good agreement among themselves). The difference between the average of all the single 
dish data (Sio)sd = (88 =t 9) Jykms"^ and that of the two interferometric measurements 
(Sio)iNT = (65±7) Jykms~^ is significant at a ~ 2a level. If confirmed at a higher significance 
level, this discrepancy may result from the presence of a colder and more extended gas phase 
that is missed in typical interferometer maps because of its low brightness (and the resulting 
low S/N per beam), rather than the lack of short baselines. The dust content of such a gas 
phase, present even in some luminous starburst/QSOs, has been detected via its sub-mm 
continuum emission (e.g. Papadopoulos & Seaquist 1999a), where dust with Tjust ~ 15 K, 
concomitant with HI and low surface brightness CO J=l-0 emission has been found. Sub- 
mm continuum imaging and multi-transition CO observations of LIRGs have shown that 
unless one has sufficient angular resolution to separate the warm, star-forming gas and dust 
from the typically, more extended non star-forming phase, the global dust continuum and 
CO line ratios will be dominated by the warm gas (e.g., Papadopoulos & Allen 2000). In 
this paper we assume that such an extended phase of non star-forming molecular gas is not 
present in Mrk231. 

3. The state of the moleculcir gas in Mrk231 

The combination of the large number of CO and HCN line detections for the ULIRG / QSO 
Mrk231 (Table 2) and the upper limit on the ^^CO/^^CO (J=2-l) line ratio of R21 > 37 
(Glenn & Hunter 2001) offers an excellent opportunity to place strong constraints on the 
state of its molecular gas, with emphasis on the dense gas phase. To do this, we used a 
Large Velocity Gradient (LVG) code, based on work by Richardson 1985, that searches a 
large grid of (n, Tk,Ax) values, where Ax = rx/(dV/dr) (rx = [X/H2], dV/dr: cloud velocity 
gradient, X: the molecule used) and locates the = Xli V^?[I^i ~ Robs.i]^ minima, where 
Robs.i and R; arc observed and model line ratios, and (J\ is the measurement uncertainty. A 
parameter space of n = (10^ — 10^)cm~^, Tk = (15 — 150) K, and an Ax range corresponding 
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to Kvir ~ 0.05 — 500 (for the standard CO and HCN abundances) was searched, where 



(e.g. Papadopoulos & Seaquist 1999b; Goldsmith 2001) indicates whether a viriahzed gas 
phase (Kvir ~ 1 within factors of 2-3), or a non- viriahzed one (Kvk ^ 1) is responsible for the 
molecular line emission (a ~ 1 — 2.5 depending on the assumed cloud density profile, Bryant 
& Scoville 1996). Strictly speaking, values of Kyi^ <^ 1 are not physical (i.e. gas motions 
cannot be slower than those dictated by self-gravity), however they can be used to indicate 
the possibility of enhanced molecular abundances instead (so that the revised Kvir can reach 
at least ~ 1). Additional constraints to the modelling are provided from estimates of the dust 
temperature - Tdust(Mrk231) ~ (47 - 54) K (Glenn & Hunter 2001; Gao & Solomon 2004a) 
- and the assumption that photoelectric and/or turbulent gas heating (and its cooling via 
atomic/molecular lines rather than continuum emission) results in Tk > Tdust (e.g. Wilson 
et al. 1982; Tielens & HoUenbach 1999). The chosen range for Tk therefore encompasses 
that expected for the molecular gas in Mrk231, with the temperature of C^-cooled, Cold 
Neutral Medium HI gas setting its upper hmit. 

The results of single-phase modeling make it clear that the CO J=4-3 and J=6-5 trace 
a different gas phase from that traced by the J— 1-0, 2-1, 3-2 transitions, i.e. we find 
no region of the (Tk, n. Acq) parameter space compatible with both the observed values of 



r2i = (2-l)/(l-0), r32 = (3-2)/(l-0), R21, anc? the high r65/43 = (6-5)/(4-3) ~ 0.66±0.26 



ratio (x^ ^2.5 over the entire parameter space). The LVG solutions derived solely from r2i, 
r32 and R21 converge to Tk - (55 - 95) K, n ~ 10^ cm'^ (best fit when [i^CO/^^CO] = 100, 
and for Tk = 75 K; = 0.7), and Kvir ~ Iba'^^^ (for rco = [CO/H2] ~ 10"^). This warm, 
non-virialised, gas phase with rio(^^CO) ~ 1 is typically found in ULIRGs, a possible result 
of strongly evolving dynamical and far-UV-intense environments (Aalto ct al. 1995), but has 
expected ratios of T43 < 0.35 and tq^ < 0.04, which are > 2 (for r43) to > 11 (for tq^) times 
lower than those observed in Mrk231. This result is illustrated in Fig. 3, where we show 
the observed CO line luminosities (in solar units) together with the values expected from 
the best-fitting LVG model derived from the lowest three CO lines and the ^^CO J=2— 1 
line only. The fluxes of the high- J lines far exceed the values predicted from this model, 
demonstrating the presence of a distinct gas component, the emission from which dominates 
the observed line fluxes. 
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3.1. The dense gas phase 

The presence of a massive gas phase in Mrk 231 that is much denser than that dominating 
the low- J CO hne emission can be inferred from its luminous HCN J=l-0 line which makes 
this galaxy stand out even amongst ULIRGs as the one with the highest HCN/CO J=l-0 
luminosity ratio (~ 0.25, Solomon et al. 1992a). The detection of HCN J=4-3 (ucrit ~ 
8.5 X 10^ cm^'^) certainly corroborates this, however a single line is not sufficient to point 
uniquely to the presence of gas with n > Ucrit- Much lower densities are still possible if sub- 
thermal excitation and/or radiative trapping (due to their considerable optical depths) were 
to be significant. Intensity ratios of widely-spaced HCN transitions are excellent probes of 
the dense star-forming gas, and provide the key to differentiating between different excitation 
processes, revealing a considerable range of properties even in starbursts with similar FIR 
and low- J CO luminosities (Jackson et al. 1995; Paglione, Jackson, & Ishizuki 1997). 

The r43(HCN) ratio measured in Mrk 231 indicates a sub-thermally excited HCN J=4-3 
line which is well below that of Arp 220 (~ 0.8; Greve et al. 2006), another archetypal ULIRG 
often used as a typical template for high-z starbursts. The physical conditions compatible 
with the dense-gas-dominated r43(HCN) and HCN(l-0)/CO(6-5) (= Rhcn/co) brightness 
temperature ratios were explored with our LVG code using HCN collisional rates for the first 
11 levels taken from the Leiden Atomic and Molecular Database LAMDaJ^I The value of 
R-HCN/co helps discriminate over the considerable range of conditions compatible solely with 
r43(HCN), while two additional constraints can be set by a) assuming that the HCN emission 
emanates from virialized gas "cells" (Kvir ~ 1; as is the case for the dense star-forming gas in 
the Galaxy), and b) by stipulating that Tk > Tjust- The best solution ranges found are those 
with Tk = 40 — 45 K and Tk = 50 — 70 K (see Table 3), though the latter satisfies Kvij. ~ 1 
only for a ~ 5 — 10 times higher [HCN/H2] abundance than that considered typical for the 
Milky Way. In starburst environments values of Tk = 50 — 70 K may be possible because of 
a higher ionization fraction of the molecular gas (Lepp & Dalgarno 1996), or higher C, C^ 
abundances (which favor HCN production, Boger & Sternberg 2005) deeper inside molecular 
clouds caused by a more vigorous turbulent diffusion of their atom-rich outer layers inwards. 

The mass of the HCN-emitting gas phase can be estimated in a manner similar to that 
used to determine the total molecular gas mass (using the ^^CO J=l-0 line), since the same 
arguments about line emission from an ensemble of self-gravitating, non-shadowing (in space 
or velocity), clouds apply. Following Gao & Solomon (2004a), 



■^http://www. strw.leidenuniv.nl/moldata/ 
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M_(H.)«2.1«(5,^)l„™. (3, 

where Tb and Lhcn are the area/velocity-averaged brightness temperature and the hne 
luminosity of an optically thick HCN J=l-0 line (the case for all solutions in Table 3). 

For the two sets of conditions that best fit the dense gas line ratios, the coefficient in the 
equation above becomes Xhcn ~ (19 — 20) (Kkms~^ pc^)~-^ (for Tk = 40 — 45 K) and 
Xhcn~ (8-9)Mo (Kkms-^pc2)-i (for Tk = 50- 70K). 

The HCN J=l-0 luminosity of Mrk231 is estimated using 

L.= / / TbdadV^ (^) f S.dV, (4) 

where AV and Ag are the total linewidth and the area of the emitting source (where in this 
case x=HCN(l-0)) respectively. Substituting and converting to astrophysically useful units 
yields 



= 3.25 X 10^ (1 + z)-^ f ^) ( ^) ' ( ^K^l K km s"^ pc^ (5) 
^ ^ VGHz/ VMpcy VJykms-V ^ 

For Dl(z=0.042)=183.4 Mpc and z/^,rest = 88.63 GHz (HCN J=l-0) the HCN J=l-0 velocity- 
integrated flux density (Table 2) yields Lhcn = 2 x 10^ Kkms~^ pc^. For the range of Xhcn 
values derived previously, this corresponds to Mdense(H2) ~ (1.6—4) x 10-'^'^ Mq, while the total 
H2 gas mass estimated from the CO J=l-0 line luminosity of Leo = 6.9 x 10^ Kkms~^ pc^ 
(^x,rcst = 115.27 GHz, and the CO J=l-0 flux in Table 2), and a standard Galactic conversion 
factor of x[fo^' = 4M0 (Kkms-i)-\ is Mtot(H2) ~ 3 x IO^^Mq. Thus, at least ~ 50% and 
maybe all of the molecular gas mass in the ULIRG/QSO Mrk231 is dense f> W^cm^^), 
quite unlike the state of the bulk of the molecular gas in our Galaxy and in low intensity 
starbursts where (n(H2)) ~ 10^ — 10^ cm"'^ (e.g. Paglione et al. 1997; Weiss, Walter, & 
Scoville 2005). 

The dynamical mass within a radius of R ~ 1.7 kpc, estimated from high-resolution 
CO images, is Mdyn ~ 3.25 x 10^° Mq (Downes & Solomon 1998, corrected for the adopted 
cosmology). Thus, either all that mass is molecular gas, or Mtot(H2) is overestimated by 
adopting an X^^^^' factor. Extensive studies of ULIRGs suggest the latter because in these 
extreme systems the molecular gas phase encompasses significant amounts of non-gaseous 
mass (i.e. stars) and thus one of the main assumptions underlying a standard Galactic Xco 
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factor (that of an ensemble of self-gravitating molecular clouds) breaks down. These studies 
find Xco ~ 1/5 X[5)'''^ (Solomon 1997; Downcs & Solomon 1998), which for Mrk231 yields 
Mtot(H2) ~ 6 X 10^ Mq, making it ~ 2.5 — 6.5 times smaller than Mdense(H2) estimated 
from Equation 3. This is clearly impossible, and argues for a revision also of the Xhcn 
values derived from Equation 3. Such a revision can be understood in much the same terms 
underlying that of the X^^'-* factor, since the dense gas phase seems responsible for both 
HCN and most of the CO J=l-0 hne emission. It must be noted that this is not done for 
many HCN-deduced dense gas masses in ULIRGs reported in the hterature (e.g. Solomon 
et al. 1992a) and thus these could be systematically overestimated by a factor of '^5. 

Interestingly, by adopting the same correction factor of ~ 1/5 for Xhcn yields Mdciisc(H2) ~ 
(3 — 8) X 10^ Mq for Mrk231 which, for its starburst- related IR luminosity of ~ 2/3 Lir 
(Downes & Solomon 1998) and Lir(8- lOOO/zm) = 3.6 x 10^^ Lq (Sanders et al. 2003), raises 
the star formation efficiency to esF = LS/Mdense(H2) ~ (300 - 800) Lq/Mq, effectively 
bracketing the maximum value of ~ 500 Lq/Mq expected from O, B, star radiation-feedback 
effects on the accreted dust /gas in star- forming GMCs (Scoville 2004). These may therefore 
be the hallmark features of the star-forming molecular gas phase in galaxies - densities of 
n > 10^ cm~^, and star formation efficiencies of ^ 500 Lq/Mq - which in Mrk231 amounts 
to most of its molecular gas mass. 



3.2. The diffuse gas phase 

Unlike the HCN J=l-0, 4-3, and CO J=6-5 line emission which are dominated by the 
dense gas, the lower-J CO lines have contributions also from a diffuse phase. The observed 
CO (6-5)/(l-0) line ratio can then be expressed as 

_ Cba (b) 
^65 — -, , ^ ^65 ' 

where Cba = fbaT^oVT^o^ expresses the contribution of emission from the dense phase (b) to 
that from a diffuse and more extended phase (a), with fba < 1 being their relative geometric 
filling factor, and T^q'^^ their CO J=l-0 brightness temperatures (all quantities are veloc- 
ity/area averages). All other CO line ratios where both phases contribute to both transitions 
can be expressed as 

(a) p (b) 

rj+ia = r ^''^ J + 1 = 1, 2, 3, 4. (7) 

J- + ^ba 
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Typically we find Cba ~ 0.6 — 0.7 which, along with the known range of Tj^^j values 
(obtained from the dense gas properties outlined in Section 3.1), allows the subtraction of 

(a) 

the dense phase contribution from the observed ratios and the estimate of the rj^^ j values. 
These and the R21 > 37 ratio (which we assume to be dominated by the diffuse phase) are 
then used as inputs into our LVG code. This in turn yields n ~ 10'^ cm^^ (Tk ~ 45 — 85 K) 
for the diffuse phase (a), with still lower densities of n ~ 300 cm~^ but Tk ~ 80 — 140 K 
also possible. In all cases Kyir > 1 (reaching as high as ~ 30), which suggests the presence 
of highly unbound gas whose large dV/dr values and high Tk are responsible for its modest 
CO J=l-0 optical depths (tio ~ 0.5 — 1.5). This phase could be confined around individual 
GMCs "enveloping" their much denser self-gravitating regions, or could be distributed over 
very different scales. Interferometric CO line imaging reveals a disk distribution with much 
of the lower-density molecular gas in the outer regions (Downes & Solomon 1998). In either 
case, the diffuse non self- gravitating gas phase is a minor contributor to the total molecular 
gas mass of Mrk231. 



3.3. Thermal balance of the molecular gas in Mrk 231 

The temperature of the molecular gas in a galaxy is determined by the equilibrium 
between heating and cooling processes. Gas heating is dominated by the photoelectric effect 
on dust grains and polycyclic aromatic hydrocarbons (PAHs), as discussed by e.g., Wolfirc et 
al. (1995). Cooling proceeds through line radiation which, in normal galaxies, is dominated 
by the [Cll] 158 /xm line (e.g., Wolfire et al. 1995, 2003; Kaufman et al. 1999), which can 
carry up to ~ 0.1 — 0.5% of the total far-infrared (FIR) luminosity of a galaxy (e.g., Stacey 
et al. 1991; Malhotra et al. 1997; Leech et al. 1999; Pierini et al. 1999; Negishi et al. 2001). 
In thermal equilibrium, the heating and cooling rates are balanced: increased heating, such 
as that resulting from enhanced star formation will be compensated by increased cooling. 
This has been used to estimate star formation rates in nearby galaxies using the observed 
[Cll] line luminosity (e.g., Crawford et al. 1985; Stacey et al. 1991; BoseUi et al. 2002). 

Measurements obtained with the Long Wavelength Spectrograph (LWS) on the Infrared 
Space Observatory (ISO) have shown that the [Cll]/FIR flux ratio decreases in galaxies of 
very high FIR luminosity, i.e., the increase in [Cll] line luminosity is no longer proportional 
to the FIR luminosity. This effect was first discovered in deep ISO measurements of a small 
sample of nearby ULIRGs (Luhman et al. 1998), where the observed [Cll]/FIR luminosity 
ratio was found to be smaller than ~ 0.05% in some of the most FIR-luminous galaxies - 
a value more than factor of 10 lower than that observed in less luminous galaxies. Various 
explanations for this effect have been proposed, including dust absorption and saturation 
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effects, self- absorption in the [C ll] line and more subtle explanations related to the detailed 
physics of dense photon dominated regions (PDRs) (Luhman et al. 2003). 

The high luminosity of the [C ll] line suggests that the hne should be observable (red- 
shifted into the submiUimeter regime) out to very high redshifts (e.g., Loeb 1993; Stark 1997; 
Suginohara et al. 1999). Indeed, ALMA will be able to detect the [C ll] emission from a Milky 

Way- type galaxy out to ^ ~ 5 (Van dcr Werf & Israel 1996). It is therefore imperative to 
understand the physical conditions determining the strength of the [C II] line more fully. 
The first attempts to detect this line in the distant Universe were carried out by Isaak et al. 
(1994) on the z = 4.7 QSO BR 1202-0725, with an upper limit implying a [Cll]/FIR ratio 
lower than in the nearby starburst galaxy M82. In a deeper integration on this same object, 
an upper limit on the [C ll] luminosity was found which is less than 0.06% of the total FIR 
luminosity (Van der Werf 1999). Recently, the first detections of [Cll] at high redshift have 
been obtained in the z = 6.42 QSO SDSS J1148+5251 (Maiolino et al. 2005), and in the 
z = 4.7 QSO BR 1202-0725 (lono et al., 2006) where the [Cll] luminosity is 0.02 - 0.04% 
of the total far-infrared luminosity. All these results are consistent with ISO observations of 
local ULIRGs, and suggest that the same physical mechanisms are responsible. 

Our present results for Mrk 231 amount to the most complete picture of the molecular 
gas attained for a local ULIRG enabling us to analyze the cooling budget of its molecular 
medium in detail, and thereby shed light on the [C II] problem. The results are summarized in 
Table 4 where we have used the mean values from our LVG models to calculate luminosities 
for CO lines which have not been observed. The HCN lines (not relevant for the cooling 
budget) resulting from the same models are presented in Table 5. The [Cll] 158 /xm line 
luminosity has been derived from the line fiux of 3.2 ±0.4- 10~^° Wcm~^ measured using the 
ISO LWS by Luhman et al. (1998). The coohng due to the 370 //m and 609 //m [Cl] fines 
was estimated from a measurement of the latter (Gerin & Philfips 2000) and by assuming 
LTE (the [Cl] energy levels are expected to be fully thermalized for the dense gas phase). 
The Tk range of the best two LVG solution ranges (Table 3) constrains the CI(2-l)/(l-0) 
brightness temperature ratio to ~ 0.8 — 1.2, for which we adopt the mean value of ~ 1 (which 
is actually measured in M82; Stutzki et al. 1997). We also list the total luminosity L^^nse 
in CO lines up to J=10— 9 for the dense phase only, as well as those the total luminosity 
Ldiffuse in CO lines from the diffuse phase, calculated from the difference between observed 
fluxes and the modeled fluxes from the dense phase. Line strengths are presented in flux 
units {Siine in Jykms"-*^), in luminosity units (Lune in Lq), and in L[i^^ luminosities (in 
Kkms~^pc^, estimated using Eq. 5); the last quantity scales directly with the intrinsic 
brightness temperature of the line, and is thus constant for thermalized optically thick lines 
originating from the same medium. 
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Table 4 shows the remarkable result that the coohng luminosity for the dense molecular 
gas in Mrk231 considering only CO approaches that of the [Cll] line, with Lco/-^fir = 
1.1 ± 0.2 ■ 10"^ while L[cn]/LFiR = 1-5 ± 0.2 ■ 10-^ for Lfir(40 - 400 /xm) ~ 2.3 ■ 10^^ ^0 
in Mrk231 (Sanders et al., 2003 for the cosmology adopted here). The situation is further 
illustrated in Fig. 4, where we show a bar histogram of the cooling lines of the interstellar 
gas in Mrk231. Error bars are observational errors for the luminosities based on measured 
fluxes, and indicate the range allowed by the LVG models for the remaining values. For the 
latter values we adopted a minimum uncertainty of at least 30% based on the measurement 
errors of the CO J— 6— 5 and 4—3 hnes. The error bars thus give a good indication of the 
expected ranges taking into account uncertainties in both the observations and the LVG 
modeling. It is instructive to compare these results to the corresponding results for the 
Milky Way, where global CO line fluxes have been measured using the COBE data, and 
converted into luminosities using a model for the spatial distribution of the emission within 
the solar circle by Wright et al. (1991). For the Milky Way the total CO cooling (relative 
to FIR luminosity) is Lco/-^fir ~ 2 • 10"^, while for [Cll] L^cui/Lfir = 2.8 ■ 10~^, i.e., the 
coohng is totally dominated by [C ll] with neghble CO cooling. The importance of the CO 
hne cooling is also reflected in the CO hne ratios: the flux ratio CO J=4— 3/J=2— 1 (which 
has been observed both objects) is 1.5 in the Milky Way (Wright et al., 1991), but 5.4 in 
Mrk231. We can also compare these results to a recent survey in mid- J CO lines and the 
two [C l] lines by Bayet et al. (2006) of a small sample of local starburst galaxies. In these 
objects the mid-J CO lines are also found to be much stronger than in the Milky Way, with 
cooling by CO lines significantly exceeding that from [Cl]. However, CO cooling in these 
objects is still insignificant compared to [Cll], in contrast to the much more extreme results 
obtained here for Mrk231. 

The cooling associated with the diffuse molecular gas (contributing only to the lowest 
three CO transitions) of Mrk231 is relatively unimportant. Thus it is natural to identify 
the dense phase with the actively star forming gas, which experiences the strongest heating, 
and must therefore also cool efficiently. Unlike the Milky Way and lower intensity starbursts 
in Mrk231 this phase contains most of the molecular gas mass, forming stars at what is 
thought to be a maximum efficiency (Section 3.1), at rates found only in dense star-forming 
cores in the Milky Way. 

The resulting picture of a dense and dominant (in terms of mass) molecular gas phase, 
emitting strongly in CO J— 6— 5 but with a suppressed [C ll] hne, points towards dense PDRs 
(e.g., Kaufman et al., 1999). In such PDRs, the high density leads to a high formation rate 
for CO (proportional to n^), while the CO photodissociation rate is less strongly enhanced 
(proportional to n). As a result, the ionized carbon layer is thin, leading to a suppressed 
[Cll] line. In addition, the transition from ionized and atomic carbon to CO now takes 
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place closer to the source of heating, resulting in a large column density of warm molecular 
gas with strong emission in the mid-J CO lines. Generalizing to ULIRGs as a population, 
our Mrk231 result suggests that the suppressed [Cll] cooling in ULIRGs is a result of high 
densities for most of their molecular gas, bathed in strong far-UV radiation fields, quite 
unlike lower intensity starburst or quiescent galaxies where much more modest amounts of 
molecular gas reside in such a high density phase. 

Our explanation predicts strong emission in CO lines up to J=10— 9 in ULlRCs. While 
such high-J transitions are not observable from the ground, they will be accesible with HIFI 
onboard the Hcrschcl satellite. Observations of high excitation CO lines from space or the 
ground, as well as careful estimates of the fraction of molecular gas mass with n > 10^ cm~^ 
in local ULIRGs using HCN transitions, will provide a critical test of the hypothesis that 
dense PDRs lie at the heart of the [C ll] cooling problem in ULIRGs. 



The wealth of molecular line and dust continuum data used to constrain the state of 
molecular gas in Mrk231 and assign most of its mass to a dense phase is rarely available for 
objects in the distant Universe. Typically, two high-J CO hues are detected (Solomon et 
al. 1992b; Solomon & Vanden Bout 2005 and references therein) , with a similar sparsity of 
dust continuum measurements. The few frequencies observed (usually at the Rayleigh- Jeans 
part of a SED) provide poor estimates of T^ust, which in turn provides little constraint on 
the range of Tk. In most cases, the often (but not always) large CO line ratios of high-J lines 
measured in high redshift objects are interpreted in terms of a warm phase, yet infact they 
trace a wider range of physical conditions. This degeneracy can be easily demonstrated by 
setting the high tq^/^s ratio measured in Mrk231 as the only constraint available. It can then 
be seen (Table 6) that besides the typical warm gas phase with moderate optical depths, gas 
as cold as Tk = 15 K but ~ 30 — 100 times more dense also reproduces the high r65/43 ratio 
that we have observed. In such a phase the large CO line optical depths, and the resulting 
radiative trapping, thermalize transitions up to high J levels so that 



(Tj_|_ij = hi/j+ij/ke), which for e.g. Tk = 17K yields r65/43 ~ 0.66 (the observed value). 

The degeneracies alluded to above reflect real conditions found in GMCs: warm, dense 
gas associated with star formation, and dense-but-cold gas in regions with no such significant 



4. Moleculeir lines at high redshifts: caveats and expectations 
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activity. In ULIRGs a massive, dense but cold phase is unlikely but not impossible during 
the rapid dynamic evolution expected for mergers (e.g., Aalto 2005). Observing at least one 
^'^CO transition is pivotal to discriminating between cold dense gas with an optically thick 
CO J=l-0 line and a less dense warmer phase where rio(^^CO) < 1. Indeed, the cold/dense 
solutions in Table 6 correspond to ^^CO/^^CO J=l-0 intensity ratios of Rio ~ 1 — 2, much 
lower than those even in the coldest GMCs (~ 5). Interestingly, setting Kyir > 1 as a 
constraint selects the conditions with Tk > 40K (Table 6), and thus may also be useful in 
"breaking" such degeneracies, though it will always depend on the assumed abundances. 

4.1. Moleculcir gas SLEDs for Mrk231: towcirds establishing local benchmcirks 

High resolution imaging of molecular lines and their relative strengths with the next gen- 
eration of mm / sub- mm arrays holds the key to unobscured views of deeply dust-enshrouded 
star forming regions in galaxies, their dynamical masses, and the molecular gas fueling the 
embedded star formation. 

This is apparent now that CO-bright H2 gas has been detected out to z ~ 6 (Walter 
et al. 2003), and distant dust-enshrouded optically faint starbursts, responsible for building 
large fractions of the stellar mass at present epochs, have been discovered (e.g. Small, Ivison, 
& Blain 1997). The need for local molecular line SLEDs for the bulk of the molecular gas in 
galaxies (and not just for their starburst sub-regions) is underlined by the fact that currently 
there are more detections of CO J=4-3 or higher-J lines in high redshift rather than local 
starbursts (e.g. Solomon & Vanden Bout 2005). 

The compact distribution of the CO, HCN emission in local ULIRGs and their large 
IR and molecular hue luminosities makes them ideal objects for estabhshing such SLEDs - 
single pointings encompass all the H2 distribution - a situation mirrored in high-z galaxies. 
The extremc-starburst/QSO present in Mrk231 makes this galaxy in particular an often- 
used benchmark for similar high-z systems (e.g. Wagg ct al. 2005). We use the constraints 
on its dense gas properties to derive the expected strengths of the unobserved CO and HCN 
transitions. The resulting CO and HCN SLEDs can then be used to assess the capabilities 
of the next generation cm (EVLA, SKA) and mm/sub-mm (ALMA) arrays, as well as the 
spaceborne Herschel Observatory in detecting star-forming molecular gas throughout the 
Universe (e.g. Van der Werf & Israel 1996). The line fluxes, calculated for a range of 
redshifts, are presented in Fig. 5, where we also show the analogous values for the Milky 
Way, with luminosites again from the COBE data (Wright et al., 1991) for all lines except 
CO J=l— 0. The luminosity of the latter line in the Milky Way was estimated using the 
CO J=2— 1 luminosity from Wright et al. (1991) and the CO 2—1/1—0 line ratio in the 
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inner 2.5° of the Milky Way as measured by COBE, which is within the errors equal to the 
global line ratio over the inner Galaxy (Fixsen, Bennett, & Mather, 1999). The Milky Way 
and Mrk231 bracket the two extremes of quiescent and starburst galaxies. A luminosity 
of LpiR = 2 ■ IO^^Lq (Wright et al., 1991) makes the Milky Way detectable with ALMA 
out z ~ 1, using CO J=3— 2 in the 211 — 275 GHz band. Higher GO lines will not be 
detectable however, though redshifted [G ll] will be observable with ALMA out to 2; ~ 4 — 8 
with significant integration time. In the case of a Mrk231-like object, the situation is quite 
different, with transitions up to GO J—9—S detectable with ALMA out to 2; ~ 10. The 
two [G l] hnes can be detected out to z ~ 2.5, and the [G ll] line will be readily detectable 
at 2; ~ 4 — 8. The high-density tracing HGN lines will typically be detectable only out to 
z ~ 1.5, where the 84 — 116 GHz band of ALMA will be the most suitable frequency range. 

It is instructive to investigate how existing detections of GO lines at high redshift fit on 
the Mrk231 template. We use the comprehensive compilation by Solomon & Vanden Bout 
(2005; their Appendix 2) of all detections of (sub)millimeter lines from high-z galaxies to 
2005, supplemented with more recent data by Riechers et al. (2006). In Fig. 6 we plot again 
the GO cooling budget of Mrk231 but now overlay the observed line luminosities of all high- 
z objects with at least two detected GO lines (one of them being GO J=3— 2 or J=2— 1). 
Ideally all fluxes should be scaled relative to one common line, but unfortunately there is 
no line which is detected in all (or even most) objects. However, almost all high- 2; objects 
with GO detections have either the J=3— 2 or the J=2— 1 hne detected (but, remarkably, so 
far never both). We therefore scale the observed line fluxes to either the GO J=3— 2 or the 
GO J=2— 1 luminosity of Mrk231 and show the results in Fig. 6. Inspection of this figure 
reveals a considerable range in J=4— 3/2— 1 line ratios, and it is clear that interpolation 
to a fictitious J=3— 2 value would have introduced significant uncertainties. Therefore we 
prefer to normalize to either the J=3— 2 or the J=2— 1 line, which, although blurring the 
comparison somewhat, has the advantage of being based on measurements rather than an 
uncertain interpolation. Physically, this procedure is acceptable since these lines have similar 
contributions from the dense gas component (70% for the 2—1 line, 80% for the 3—2 hne) 
in our Mrk231 fiducial model. 

Figure 6 shows that Mrk231 is a reasonable template for most high- 2; objects that have 
been detected in GO, with line ratios within a factor of two from those of Mrk 231 up to about 
J=6— 5. This does not imply that all high- 2; objects would have similar line ratios: there 
is clearly a selection effect in the sense that high-2; objects are followed up with low-J GO 
observations only after having being detected via their luminous high- J line emission. With 
this caveat in mind, Mrk 231 appears to be a useful template for high- 2; galaxies detected in 
GO to-date. It is interesting to note that this is the case also for its dust continuum emission 
(see e.g.. Fig. 2 of Blain et al., 2002). 
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Two notable outliers in Fig. 6 are HRIO and APM 08279+5255. HRIO is an extremely 
red object at a modest redshift of z = 1.439 which, while a ULIRG (LpiR ~ 6 ■ 10^^ Lq), 
has line ratios indicative of low excitation conditions more similar to the Milky Way than 
to Mrk231 (Papadopoulos & Ivison, 2002). Detections and imaging of such low excitation 
objects beyond 2; ~ 4 (even with such large Lpin) will be difficult even with ALMA. At the 
other extreme, APM 08279+5255 shows unusually high excitation, with CO line luminosities 
rising all the way up to at least J=10— 9: gas at high gas densities is indicated by the 
detection of luminous HCN J=5— 4 by Wagg et al. (2005). While these authors can fit their 
data with a single gas phase of n(H2) ~ 4 x 10^ cm~^, this object is strongly lensed (Lewis et 
al., 2002), and so it is hkely that the observed fluxes are dominated by a selectively amphfled 
region of high-density molecular gas. Another indication of a singularly different molecular 
line excitation of this distant starburst/QSO is its HCN(5-4)/CO(4-3) luminosity ratio of 
L;i„^(HCN)/L;i„^(CO) ~ 0.21 - 0.36 (from Wagg et al. 2005 and Downes et al. 1999), while 
from the observed CO J=4-3 and the expected HCN J=5-4 luminosity (deduced from the 
best LVG solutions for the dense gas) the much lower values of ~ 0.045 — 0.066 are derived 
for Mrk231. 

Agreement with (or deviation from) the Mrk231 SLED appears to depend primar- 
ily on the mean gas density of the phase dominating the emission. This is suggested by 
further inspection of Fig. 6 in relation to the HCN line luminosities where available. In 
IRAS F10214+4724, VCV J140955.5+562827, and the Cloverleaf quasar, HCN J=l-0 has 
been detected, but with an HCN(J=l-0)/CO(J=3-2) ratio much lower than in Mrk231, 
and these galaxies reveal somewhat lower CO line excitation than Mrk 231 as well. In con- 
trast, the HCN(J=2-l)/CO(J=2-l) ratio in SDSS J1148+5251 is higher than in Mrk 231, 
and this object shows also higher excitation in the CO lines. Construction of molecular gas 
SLEDs out to the high excitation regimes associated with the star-forming gas for a large 
number of local LlRCs will allow a more thorough examination of such differences, while 
establishing much-needed benchmarks for the interpretation of the still sparse data at high 
redshifts. 



5. Conclusions 

We report on the detections of the CO J=6-5, J=4-3 transitions of the archetypal 
nearby ULIRC/QSO Mrk 231. These first detections of high-excitation CO lines in a local 
extreme starburst, along with hereby reported CO J=3-2, HCN J=4-3 detections and CO 
J=l-0, 2-1, ^^CO J=2-l, HCN J=l-0 literature data, allow robust constraints to be placed 
on its molecular gas properties. Our conclusions can be summarized as follows: 
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1. The high excitation CO J=4-3 and J=6-5 hnes trace denser (> lO^cm"^), gas phase 
different from that dominating the lower three CO transitions. This phase is also 
responsible for the observed HCN line emission, and contains > 50% of the total 
molecular gas mass in this galaxy, quite unlike quiescent or lower intensity starbursts 
where only a few percent of the total molecular gas mass reside in such a phase. 
This dominant dense gas phase fuels the starburst in Mrk231 at almost the maximum 
expected star forming efficiency of L|p'/M(n > lO^cm"^) ~ SOOLq/Mq. 

2. This preeminence of dense and relatively warm (Tk ~ 40—70 K) gas in Mrk231 presents 
a totally different thermal balance from that found in more quiescent galaxies, with 
CO comparable to [Cll] emission line cooling. If confirmed, through observations of 
other ULIRGs, this may provide the explanation of the long-standing issue of their 
very low [Cll]/IR luminosity ratios namely, bright high-excitation CO lines from dense 
far-UV photon-dominated molecular gas are the major and even dominant coolants 
with respect to their faint [C ll] line. 

3. Since the excitation conditions probed by the mid- J CO lines and low-excitation atomic 
and ionic fine-structure lines are very similar, mid-J CO lines provide an excellent 
diagnostic complementary to e.g., PACS on the ESA Herschel satellite, and proposed 
future missions such as the Japanese SPICA project. With ALMA these diagnostics 
will be available at unprecedented spatial resolution. 

4. The well-constrained dense molecular gas SLED of Mrk231 provides a very valuable 
template for direct comparisons with starburst/QSO systems at high redshifts, with 
the caveat that it may not be representative of the ULIRG population as a whole (e.g. 
Mrk231, unlike Arp220, has a sub-thermal HCN(4-3)/(l-0) line ratio). Nevertheless 
it is a template unaffected by the strong differential lensing that may be skewing the 
intrinsic SLEDs of several high redshift galaxies towards those of a more compact, 
highly-excited, star-forming and even AGN-related molecular gas phase. 
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Fig. 1.— Mrk231: a = 12^ 56'" 14. 18% 5 = +56° 52' 25.8" (J2000). The CO J=3-2 (thick 
hne), 10x[HCN J=4-3] (thin hne) spectra at resolution Ai^^h = 25 MHz (~ 21.6 kms"^), 
with thermal rms errors of ST^ ~ 9mK (CO J=3-2), and ST^ ~ ImK (HON J=4-3). 
Velocities are relative to vi^sr — 12650 km s~^. 
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Fig. 2. — High excitation CO transitions: CO J=6-5 (thick hne) and CO J=4-3 (thin hne), 
at a common rcsohition of ~ 20kms~^. Thermal rms errors: 5T\[A — 3) ~ 12 mK, and 
(5T^(6 — 5) ~ 15 mK. Velocities are relative to i'lsr — 12650 kms~^. 
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Fig. 3. — Observed luminosities of CO J=l— 0, 2—1, 3—2, 4—3 and 6—5 (and interpolated 
value for CO J=5— 4) shown by the drawn boxes with error bars; the dashed lines give 
the expected luminosities based on the best-fitting single component LVG model derived 
considering only the three lowest CO lines and ^^CO J=2— 1. 
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Fig. 4. — Cooling budget of the molecular gas of Mrk231. Plotted are the total CO hue 
and [C ll] hne luminosities from Table 4, which represent the sum of the diffuse and dense 
gas phases. The bar for [C l] represents the sum of the expected luminosities of the two [C l] 
hnes (see text). Error bars are observational errors for the luminosities based on measured 
fluxes, and indicate the range allowed by the LVG models for the remaining values. For the 
CO lines deduced from these models we use a minimum uncertainty of 30% (based on the 
measurement errors of the CO J=6— 5 and 4—3 lines). 
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Fig. 5. — Predicted fluxes of Mrk 231 and the Milky Way as a function of redshift. The dif- 
ferent panels show the lines redshifted into 4 different frequency intervals, as indicated at the 
top of each panel. Also indicated are instrumental sensitivity limits (5 a point source limits in 
1 hour at a velocity resolution of 300kms~^; it is assumed that the full line flux is contained 
in this velocity interval). These limits have been calculated using the on-line sensitivity 
calculators for ALMA (|http: / /www.eso.org/ projects/alma/ science/bin/sensitivity.html ) us- 
ing 50 antennas, and GBT (|http : / / www. gb . nr ao . edu / GB T/ setups / senscalc . ht ml ) , and are 
projected values for the EVLA. These figures assume a fiat A-dominated cosmology with 
Hq = 71kms-i Mpc"^ and = 0.27. 
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Fig. 5. — (continued) 
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Fig. 6. — CO line luminosities in various high-redshift objects compared to the Mrk231 
template. Lines luminosities are normalized to the Mrk231 values at either CO J=3— 2 or 
CO J— 2—1. Observed hne fluxes are taken from the compilation by Solomon & Vanden 
Bout (2005), supplemented with more recent data from Riechers et al. (2006). 
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Table 1. CO J=l-0 fluxes measured for Mrk231 



Telescope 


Sline (Jykms ^) 


Reference 


IRAM 30-meter 


100 ± 20 


Solomon et al. 1997 


NRAO 12-meter 


100 ± 20 


Papadopoulos Sz Seaquist 1998 


FCRAO 


82 ±15 


Young et al. 1995 


IRAM 30-meter 


72 ±15 


Kriigel et al. 1990 


OVRO 


62 ± 10 


Bryant & Scoville 1996 


IRAM PdBI 


68 ±10 


Downes & Solomon 1998 
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Table 2. CO and HCN line fluxes and ratios 



Transition 


Sline 


Line ratio*^ 


References'^ 




(Jykms^"'^) 






CO J=l-0 


88 ±9*^ 




Table 1 (single dish data) 


CO J=2-l 


315 ± 30"^ 


0.90 ±0.13 


1, 2,3 


CO J=3-2 


568 ± 801^ 


0.71 ±0.12 


4, this work 


CO J=4-3 


980 ± 230 


0.70 ±0.18 


this work 


CO J=6-5 


1465 ± 500 


0.46 ±0.16 


this work 


HCN J=l-0 


15 ±3 




5 


HCN J=4-3 


65± 13 


0.27 ±0.08 


this work 



Velocity/area-averaged brightness temperature line 
ratios rj+i J = (Tb(J ± 1, J))/(Tb(l, 0)). 

^Averages of the values extracted from the data in the 
listed references. 

'^l. Papadopoulos & Scaquist 1998; 2. Downcs & Solomon 1998 (two 
values reported); 3. Glenn, & Hunter 2001 (CSO for S/T ~ 
50 Jy/K); 4. Lisenfeld et al. 1996; 5. Solomon, Downes, & Radford 
1992a 
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Table 3. Physical conditions of the dense gas phase 





n(H2)'^ 


Ahcn'^ 


Kvi, 


b 

r 


Rhcn/co (r43(HCN))^ 


(K) 


(cm^"^) 


(kms^^ pc^^)--*" 








15 


3 X 10^ 


3 X 10"^° 


5.9a- 


1/2 


2.67 (0.26) 


20-25 


3 X 10^ 


3 X 10-^ 


0.18a- 


-1/2 


1.97-1.55 (0.24-0.28) 


30 


3 X 10^ 


10-10 


17.8a- 


-1/2 


0.90 (0.27) 


35 


10^ 


10-9 


3.08a- 


-1/2 


0.80 (0.28) 


40-45^^ 


3 X 10^ 




0.56a 


-1/2 


0.72-0.64 (0.28-0.30) 


50-70^ 


10^ 


10-"^ 


O.la- 


1/2 


0.65-0.47 (0.23-0.28) 


75-90 


10^ 


3 X 10-^° 


10.25Q; 


-1/2 


0.25-0.21 (0.25-0.27) 


95-105 


3 X 10^ 


3 X 10-9 


1.87q;- 


-1/2 


0.21-0.19 (0.27-0.28) 


110-150 


10^ 


3 X 10-^ 


0.32a- 


-1/2 


0.19-0.15 (0.26-0.30) 



^Parameters corresponding to best LVG solutions for r43(HCN) ratio 
(x^ < 0.3; see section 3) 

''Prom Equation 2, rncN = [HCN/Ha] = 2 x 10"^ (Irvine, Goldsmith & 
Hjalmarson 1987; Lahuis & van Dishoeck 2000) (a ~ 1 — 2.5 see text). 

^The (dense gas)-dominated HCN(l-0)/CO(6-5) and HCN(4-3)/(l-0) 
brightness temperature line ratios for each LVG set of parameters. 
Observed ratios: Rhcn/co = 0.62 ± 0.24 and r43(HCN) = 0.27 ± 0.08). 

'^The best solution ranges are indicated in boldface. 
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Table 4. Cooling budget of the molecular gas in Mrk231 



Transition 


'S'line 


Tl 

-^line 


Lline 


Notes 




[Jy kms""-*-] 


[Kkms-ipc^] 


[Le] 




CO J = 1 - 


88 


6.9 ■ 10^ 


3.4-105 




CO J = 2 - 1 


315 


6.2 ■ 10^ 


2.4 ■ 10^ 




CO J = 3 - 2 


568 


5.0 • 10^ 


6.6 ■ 10^ 




CO J = 4 - 3 


980 


4.8 • lO'' 


1.5 • 10^ 




CO J — — 4 




4.0 • 10 


2.5 • 10' 




CO J = 6 - 5 


1465 


3.2 • 109 


3.4 • 10^ 




CO J = 7 - 6 




2.8 • 10^ 


4.7-10^ 


b 


CO J = 8 - 7 




2.2 ■ 10^ 


5.6 • 10^ 


b 


CO J = 9-8 




1.4- 10^ 


4.9 ■ 10^ 


b 


CO J = 10-9 




5.5 • 10« 


2.7-10^ 


b 


CO dense phase 






2.5 • 10^ 


c 


CO diffuse phase 






8.9 • 10^ 


d 


[Cl] 609 //m 


200 


8.7 • 10« 


3.3 • 10^ 




[C l] 390 fj,m 




8.7 • 10« 


1.5 • 10^ 


e 


[Cll] 158 //m 






3.4 • 10^ 





^ line strength estimated by interpolation (CO J+1 — > J, J+1 > 4 
lines are dominated by the dense gas phase). 

^ line strength estimated using the mean value resulting from the 
LVG models 

*^sum of CO line strengths from the dense phase using the mean 
of our LVG models 

'^sum of observed CO hne strengths minus modeled CO line 
strengths from the dense phase 

^ line strength estimated from the 609 /xm line (see text) 
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Table 5. Observed and modeled HCN line luminosities of Mrk231 



Transition 




'S'line 


T' 

-^line 


Lline 


Notes 






[Jy kms~^] 


[K kms~^ pc^] 


[Le] 




HCN J = 1 - 





15 


2.0 • 10^ 


4.5 • 10^ 




HCN J = 2 - 


1 




1.6 ■ 10^ 


2.9 ■ 10^ 


a 


HCN J = 3 - 


2 




9.7- 10*^ 


5.8-105 


a 


HCN J = 4 - 


3 


65 


5.5 ■ 10^ 


7.7-105 




HCN J = 5 - 


4 




2.7- 10^ 


7.4 - 105 


a 


HCN J = 6 - 


5 




9.2 • 10^ 


4.4 - 105 


a 


HCN J = 7 - 


6 




2.5 • 10^ 


1.9-105 


a 



^ line strength estimated using the mean value resulting from 
the two best ranges of LVG solutions (Table 3) 
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Table 6. Conditions compatible with the CO (6-5) / (4-3) ratio^ 



Tk 


n(H2) 


Ago 




Kvir*^ {no' ) 


(K) 


(cm^^) 


(kms^^ pc^ 






15 


3 X 10^ 


3 X 10- 


-4 


0.03a-V2 ^) 


20 


10^ 


3 X 10" 


-4 


0.16a-i/2 1) 


25 


10^ 


10-^ 




0.48a'i/2 1) 


30-35 


3 X 10^ 


3 X 10" 


-4 


0.3a-i/2 ^) 


40-45 


3 X 10^ 


10^6 




28a-^/2 (i.2_i.5) 


50-75 


3 X 10=^ 


10-^ 




0.89q;-^/2 (5_iq) 


80-135 


10^ 


3 X 10" 


-6 


16q;-V2 (0.08-0.45) 


140-150 


3 X 10=^ 


3 X 10- 


-5 


2.96q;-V2 (0.77-0.88) 



^Observed value: 105/43 = 0.66 ± 0.26 

''From Equation 2, and rco = [CO/H2] = 10"! 

^The optical depth of the CO J=l-0 line. 



